Transport coefficients in non-quasiparticle systems 

Antal Jakovac^ 

Insitute of Physics, BME Technical University, H-1 1 1 1 Budapest, Hungary 



Abstract. Transport coefficeints, in particular the shear viscosity to entropy density ratio is studied in systems 
where the small- width quasiparticle assumption is not valid. It is found that 77/ i- has no unversal lower bound, the 
minimal value depends on the system and the temperature, and can be even zero. We construct models where the 
l/4n conjectured bound is violated. 



1 Introduction 

Recent results from RHIC [1] suggest that hadronic mat- 
ter near Tc behaves as an almost ideal fluid [2]. Matter 
with small viscosity is very hard to described theoretically, 
based on model calculations. The main reason is that the 
usual perturbative approach expands the theory around the 
free gas limit, where, by construction, the elementary ex- 
citations are in the ballistic regime with infinite mean free 
path. Since all transport coefficients are proportional to 
the mean free path, we obtain infinite results for them in 
the unperturbed regime. With formulae we can state that 
77 ~ l/cr ~ l/g'^ X logs in QCD, where cr is the cross 
section, g is the coupling constant. In fact transport de- 
scription assumes diffusion, the transport coefficients are 
diff'usion constants. If we enforce a diffusive description to 
the ballistic regime, we obtain infinite diffiision constant. 

Therefore the small shear viscosity coefficient mea- 
sured at RHIC means that the elementary excitations of the 
matter after the collision have very short mean free path. 
Short mean free path means short lifetime, large width, 
and so a not a small-width quasiparticle system. If we still 
want to enforce an interacting free gas picture, it can be 
done only with very strong interactions. To describe such 
a system requires new approches. 

A possible exact method would be using lattice Monte 
Carlo (MC) simulations. One can measure {Ti2(x)Ti2(Q)) 
correlator on lattice in Euclidean time, and use exact in- 
tegral relations [3] to determine the commutator C(x) = 
{[Tiiix), 7"i2(0)]). One has to invert this relation based on 
only a discrete set of time arguments. To make it fea- 
sible one has to use prior knowledge on the result (eg. 
C{ll> > 0) > 0), and incorporate it via the Maximal En- 
tropy or similar method. Then one can extract the shear 
viscosity from the Kubo formula as 77 = lim xhe 

hard point of this analysis is that the Euclidean correlation 
function depends only very weakly on the derivative of C 
at zero frequency, ie. viscosity [4]. This means that one has 
to have large systematical uncertainties. Nevertheless one 
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can give estimates as 77/i = 0.102(56) at T = l.lATc [3], 
where the uncertainty reflects the statistical error. 

Another method to treat strongly interacting models to 
find a dual description with small coupling. For QCD such 
dual theory is not known, the "closest" theory where a dual 
theory can be constructed is the N - A- super Yang-Mills 
(SYM) theory. If the coupling and the number of colors 
are infinite {Nc » 1 and A - g^Nc » 1), then the dual 
model is a five dimensional gravity with AdS spacetime. 
Then, from graviton scattering, one can obtain the shear 
viscosity coefficient and rjjs, the latter turns out to be 1 1 An 
[5]. For a wide class of theories this is a lower bound [6]. 

Although strictly speaking there are no exact perturba- 
tive methods to treat the strongly interacting matter, one 
can give arguments to estimate 77/5 [7]. This argumenta- 
tion uses the quasiparticle picture: there rj ~ er and s ~ n 
where e is the energy density, r is the quasiparticle lifetime 
and n is the particle density. In this case rj/s ~ Et where 
E is the particle energy. For quasiparticles E > AE, and 
so rj/s > AEt > ti. Since the so obtained value is close to 
the one given by the N - A SYM at large couplings, it is 
tempting to interpret l/47r as a universal lower bound for 
any existing matter. 

This strong statement has been discussed a lot in the 
recent literature, and there are counterarguments given 
against its generic validity. In the 5D gravity models, 
one can construct higher curvature and dilaton models [8] 
where the l/47r is no more a lower bound. An explicit 
couterexample was constructed in [9], where the authors 
use a multicomponent model where the mixing entropy can 
tune the 77/i to zero. 

As the argumentation of [7] is concerned, it is strongly 
based on quasiparticle picture. In an interacting model the 
spectral function always contains other structures than the 
quasiparticle peak, and these can change the conclusions 
based on quasiparticles. For example the presence of a 
multiparticle continuum decreases considerably the 77/i ra- 
tio [10]. 

From the experimental side, one can analize the 
anisotropy of the created particles, in particular the V2 co- 
efficient. From hydro fits one can readily give an upper 
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bound for the rjls ratio [11] {rjls < 0.16), but it is much 
more difficult to present a lower bound. Taking into ac- 
count the quadratic pj dependence of t]Is [12], the true 
value can be very close to the conjectured lower bound 
[13]. 

In other, more controllable experiments it is found [14] 
that supercritical fluids may have very small rjls (or phys- 
ically equivalent) ratio. 

All these spectulations reflect the fact that the status of 
the minimum of the rjjs ratio is still not clarified. The goal 
of the present work is to determine, what can be said about 
this ratio in a generic model, without assuming a small- 
width quasiparticle picture. The generic model will be de- 
fined by giving the energy levels of the system (density of 
states, DoS). For a more detailed discussion see Ref.s [15] 
and [16]. 

The paper is organized as follows. In Section 2 we 
analyze the rjl s ratio in a generic model and determine 
the minimum of it. In Section 3 we consider some mod- 
els where the conjectured IjAn bound can be violated. We 
close the discussion with conclusions in Section 4. 



2 Transport and entropy in generic modeis 

A definite difference between the quasiparticle DoS and a 
generic one is the presence of continuum, ie. there is no 
dispersion relation, at fixed spatial momentum the system 
has still a lot of energy levels. To give an account for this 
property we define 

y InXnl = y y f-^ QQip) \p, Q) {p, Q\ (1) 

where g denotes the DoS also called spectral function, 
Q denotes conserved quantities (quantum channel), p - 
{po, p) is the total energy-momentum of the state. V is the 
volume, and we will use finite volume normalization for 
the states. 

In order to be able to study the transport coefficients in 
general, we define ffie correlators C for generic conserved 
current /, and the transport coefficient D by the Kubo for- 
mula (linear response ffieory): 

dx) = {[Jiix), JiiO)]} D = ]im (2) 

<j-»0 0) 

for the shear viscosity D - rj and 7,- — > Tn- The expecta- 
tion value is taken in equiUbrium: 

(A) = 1 Tr e-^^A, Z = Tr e"^" = e"^^, (3) 

where F is the free energy. 

To find out a generic formula for the transport coeffi- 
cients we insert an energy-momentum eigenbasis into ffie 
formula, and use ffie fact that translation is generated by 
energy-momentum operator as 7, (x) = e'^^y,(0)e~'^^. Then 



we find 

C(x) =^Yj {n\e-^"[JM,Jm\n) = 
= ^ 2((nk-^«y,(:c)|m) <m|7,(0)]|n) - 

n,m 

(mk-^^y,(0)|n)<n|/,(x)]|m)) = 
^ 1 2 ,KPn-P.)^ _ ,-fiE„^ I (nUmm) t (4) 

n,m 

By inserting the DoS we find 

X (g-^*" - e-l^"") I (k, "KUimq, Q) f-, (5) 

After Fourier transformation, wiffi p = zero spatial 
momentum 

dco, p = 0) = — 2, e-Tf (%a(k, ko + oj)x 

X [e-l^^ - I (k, ko, '?C|7,(0)|k, ko + oj, Q) p. (6) 

We can take into account ffiat ffie current cannot change 
the quantum numbers without changing the energy or mo- 
mentum, then the diffusion constant reads 

(7) 

In ffie free theory the expectation value of ffie current can 
be written in volume normalization as {k\Ji\q) = qki/(koV), 
where q is the charge carried by the current; for the shear 
viscosity q kj. In general ffie current matrix element 
is proportional to the velocity y, = ki/ko, and in volume 
normaUzation it is inversely proportional to the volume. 
So ffie generic formula is proportional to the free case, and 
we can write 

\{k,'K\Jimk,'K}f = J^(k^)(0^ , (8) 

where J'%{1(^) is a nonperturbative correction factor. Then 

As a final step we can average over the spatial angular de- 
pendence. This gives k^ - k^/3 and fc^A:^ = (k-^)'^/15. Then 
the shear viscosity reads 

After having defined the shear viscosity, we can try 
to give a formula for ffie entropy density. This is, how- 
ever, a conceptually more difficult task. The point is that 
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entropy density, being a statistical concept, is not sensi- 
ble in the microscopic theory. Indeed, if we consider the 
volume dependence of the free energy defined from the 
partition function in eq. (3), at small volumes it can show 
any volume dependence (usually grows very fast), only 
at large volumes will it be proportional to the volume. 
There is a "crossover" size L, where statistical treatment 
starts to be sensible. This size corresponds to the coarse 
graining scale, beyond that two neighboring volume ele- 
ments (grains) of linear size L interact dominantly through 
the surface, which modify the total energy of them only 
weakly. This means that interactions are effectively cut off 
at scale L, while statistics start to be valid above L. The 
value of L can be read off the linear size of the cross sec- 
tion: for a strongly interacting theory L is large, a weakly 
interacting theory has small L. 

Having said that we choose a volume V - Lr' Xo define 
the free energy density and entropy density as 



T 

f = — 7 In 



1 +L- 



V 



Q-Kik) e 



-Mo 



s = - 



So finally we obtain 



d T 
dfJ? 



dT' 
(11) 



(12) 



.3 r d'^k 



After having found this formula we can ask whether 
there is a lower bound in this formula. The generic analysis 
can be found in 1 16|, but the generic form can be easily 
understood. For that consider the small entropy case: then 
in the log in the denominator the 1 will be dominant, and 
we can write 



P y fjT^ (k^ 



small s 



From here we see that the numerator is proportional to 
while the denominator by g. If g is zero almost everywhere 
and very large at some point - this is the small width quasi- 
particle limit -, g^ is even larger at these point, and the 
result will be large. Therefore in small width quasiparti- 
cle limit rjl s is large. On the other hand, is g is small ev- 
erywhere, then g^ is even smaller, and the result is small. 
This means that small rjl s can be expected in system where 
there are no long-lived quasiparticles! 

To be more quantitiative we rewrite our formula as: 



m) 



smaUs «£>» 

where the rescaled DoS reads as 



gqcik) = 



(k^) 



2^2 



l5Zklil+pko) 



J^ik)g<Kik), 



(14) 



(15) 



and the averaging is interpreted as {{...)) = N ^ J dA{. . . ) 
where the integration measure is 



dA 



■K 



^^^'«^^^^^°^W(.).-\ (16) 



2\2 



and the normahzation factor reads 



N = \5ZT 



V r d^k 



(k2)2 



j:^\k)e-^'^. (17) 



Now we can apply the Schwarz inequaUty: {{o^)) 



{{g}} , and we see: 



smalls 



(18) 



What we see is that, although rj/s has a minimum, but it 
is no universal (like 1 /4n), it is model- and environment- 
dependent. The minimal value can even be zero by tuning 
the system to reach yV — > oo (cf. [9]), or go to i — > ie. to 
zero temperature. 

A more detailed analysis, including also systems with 
large entropy density shows [16] that 



mmimum 



T{L\s) 
Nq{LTY' 



(19) 



where T{x) ~ x for small x and ~ e^/x for large x, and 
Nq is the effective number of quantum channels (particle 
species). 



3 Model calculations 

After having found the minimal value, we can try to con- 
struct physically relevant models where this minimum can 
be reached. To this end we make some simphfications in 
the above calculation. First, we use a generalized quasipar- 
ticle systems, and substitute the DoS by the quasiparticle 
spectral function. Then the free energy reads for bosonic 
or fermionic systems 



^'^l f = Tf^,gQPik)iT)^l±e-^^). 

O O -J '^^"z 



(20) 



A second simplifiacation is that we omit the radiative cor- 
rections from the current matrix elements, ie. SF^ik) = 1. 
Third, we take only one quantum channel. In this way the 
"reduced" viscosity coefficient reads 



/3 r d*) 
15 J '(2^, 



2\2 



(k^ 



-Mo r?- 



QQpik). 



(21) 



3.1 Small width case 



To check the validity of our approximations, we apply it to 
a small-width quasiparticle case. Assume that the lowest 
lying states can be approximated with Breit-Wigner form 



QQpik) = 



ir 



{ko-BkY + r^' 



(22) 
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where is the dispersion relation. 

In the small width limit £>(?)^ ~ t 2^<5(?o - Sg). The 
formulae (20) and (21) can be calculated [16] with the re- 
sult 



- = -f(-), 



(23) 



where the function / depends on the dispersion relation, 
and m is the mass scale. For example for Sk - k, f - 
SAO/n'^ for bosons, / = 4320/(7;:'') for fermions; if = 
m+^, f = 30nT/m. 

If we consider a massless theory (eg. conformal field 
theory) then the only scale is T, and so F ~ T. This means 
that fj/s =constant. The lower limit of this constant may 
come from infinte coupling, 1 /4n. 

In the massive case the width of the particle should 
come from scattering process, therefore at low tempera- 
ture we expect F ~ e~'^^^, where M is the energy of the 

lowest scattering state. In this case Fj/s ~ Te'^^^ oo. 

This means that in the small width quasiparticle case 
there is a lower bound, which comes from the massless 
theory. 



3.2 Broad spectral function 

For an opposite case consider a flat spectral function: 



gQp(ko, k) 



2n 



£2 -£] 



■6>(£i < A:o < £2) (24) 



which is a step function, where £'i,2(fc) = ■y^k-'- + m\ . 
small temperatures (T < m\) we find [16] 



.At 



s OT2 — yn\ 



(25) 



This formula suggests that by broadening the energy distri- 
bution, the viscosity to entropy density ratio has no lower 
bound, it can be made vanish. This is exactly the message 
of the generic analysis: since q is normalized, broadening 
means small values, and so q everywhere. This re- 
sults in the smallness of r]js ratio. 



3.3 System with zero mass excitations 

Although the above example may be not physical, there 
are examples where the spectral function has no isolated 
quasiparticle peak, but it starts with a continuum. This is 
the case when an interacting system contains zero mass 
particles: then the multiparticle cut start directly from the 
quasiparticle peak. The interpretation is that a charged ob- 
ject is always surrounded by soft gauge bosons. For exam- 
ple in QED at one loop level the electron spectral func- 
tion [17], up to log corrections, is proportional to QQp{k) ~ 
(k^ - w^) ' , ie. it is divergent at the mass shell, and fall as 
l/(cL>-Sk). By BIoch-Nordsieck resummation it changes to 
l/iio - Sk)^^'', where b is the beta function. 

In order q to be normaUzable, the singular behaviour 
at the mass shell must be smoothed out by some physical 



process (cf. also [18]). Near the threshold (which is the 
dominant regime at low temperatures) we approximate 



■2\w 



QQp{q) = Cqo0iq - M)iq - M') 



(26) 



where w > -I for normalizability, and C is a constant, its 
dimension is [C] = [£']-2(i+«;) sj„^.g ,j ^ and f ~ C, 
therefore C remains in the ratio. In the massive and mass- 
less case we find 



and CT 



l{l+w) 



r-»o 



0. (27) 



In the limiting w - -\ case the rjls is constant, this is the 
case in some conformal models. 

This analysis suggests that in a system with zero mass 
exciations, but normalizable spectral functions the rjjs ra- 
tio should be vanishing at very small temperatures. 



4 Conclusions 

In this paper we discussed the behaviour of the the shear 
viscosity to entropy ratio in systems where the small-width 
quasiparticle approximation is not necessarily true. We 
have found that, although rjl s has a lower bound in each 
given system at a given temperature, but this bound is 
system- and temperature dependent; for small entropy case 
we found for the minimum ~ ,, ''4 where L is the inter- 

action length and Nq is the number of effective quantum 
channels (particle species). This minimum can be smaller 
than l/4;r. This fact was demonstrated by constructing 
models with this property: if the spectral function does not 
go to an isolated quasiparticle peak at zero temperature, 
which is a natural setup in systems with zero mass excita- 
tions, then 77/i — > at zero temperature. 
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